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1.0  SUMMARY 

Predicting  the  performance  of  a sonar  system  that 
incorporates  lofar  processing  requires  quantitative  infor- 
mation about  the  detection  performance  of  the  lofar.  If 
such  information  is  not  immediately  available  for  a particular 
model,  it  would  be  useful  to  have  a means  for  predicting 
its  performance  given  the  processing  and  display  parameters. 

A means  for  doing  this  is  developed  in  this  report.  The 
results  would  also  serve  as  a display  design  tool  and  as  a 
means  for  determining  ways  to  utilize  equipment  more  effectively. 

Lofar  equipment  permits  a sonar  operator  to  detect  a 
narrowband  signal  in  a noise  background.  The  equipment  consists 
of  a frequency  analyzer  and  a two-dimensional  display  that  is 
intensity-modulated  (or  reflectance-modulated  in  the  case  of 
electrosensitive  paper  recorders).  In  essence,  the  display 
presents  a time  history  of  the  frequency  distribution  of  the 
input  energy.  A steady  narrowband  signal  appears  as  a line  of 
different  intensity  from  the  background. 

The  performance  that  can  be  achieved  with  lofar  equipment 
depends  heavily  on  the  performance  of  the  human  visual  process- 
ing system  and  on  how  well  the  display  format  is  matched  to  it. 
This  system  has  been  the  subject  of  study  for  many  decades  by 
many  researchers  of  diverse  disciplines.  The  body  of  research 
found  to  be  relevant  to  this  problem  pertains  to  the  subject 
of  receptive  fields.  Some  of  this  research  is  reviewed  in 
Section  2.0  of  this  report.  In  addition,  psychophysical  data 
were  found  and  employed  to  characterize  the  type  of  receptive 
field  whose  structure  is  appropriate  for  the  detection  of 
elongated  signal  traces. 
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Based  on  the  characteristics  of  that  type  of  receptive 

field,  the  detection  performance  was  derived  by  the  methods  of 
statistical  detection  theory.  The  analysis  (Section  3.0) 
consists  of  four  major  parts: 

1.  Derivation  of  selected  statistical  measures  of 
the  lofar  display  markings. 

2.  Description  of  the  retinal  image  of  the  display 
surface  pattern. 

3.  Derivation  of  selected  statistical  measures  of  the 
response  of  the  receptive  field  to  the  light  pro- 
ducing the  retinal  image. 

4.  Calculation  of  the  detection  performance  from  the 
derived  response  measures. 

The  major  result  of  the  derivation  is  a pair  of  simple 

formulas:  One  (80)  for  the  detection  parameter  from  which  the 

probability  of  detection  can  be  calculated;  and  another  (8l)  for 
the  signal-to-noise  ratio  required  for  a specified  probability 
of  detection.  (Ten  times  the  logarithm  (to  the  base  10]  of 
this  ratio  is  the  recognition  differential.) 

The  results,  which  apply  to  a fully-alerted  operator, 
depend  on  two  sets  of  parameters: 

1.  Display  Parameters 

• Length  and  width  of  the  display  surface 

• Number  of  frequency  cells 

• Number  of  statistically  independent  time  samples 


2 


Report  No.  3225 


Bolt  Beranek  and  Newman  Inc. 


2.  Viewing  Parameters 

• Angle  and  distance 

Application  of  the  results  to  an  example  with  representative 
parameter  values  predicts  the  following  performance  characteristics 

1.  Performance  improves  with  record  length,  but  at 
a diminishing  rate. 

2.  Performance  improves  if  the  viewing  angle  is  decreased. 

3.  Performance  improves  with  the  width  of  a frequency 
cell  on  the  display,  but  only  to  the  point  where 

the  cell  width  matches  the  width  of  the  receptive  field. 

4.  A modest  advantage  can  be  obtained  by  compressing 
the  display  in  the  time  dimension.  The  advantage 
obtained  depends  on  signal  duration  and  viewing  angle. 

Although  the  equations  apply  to  a single  look  or  eye 
fixation,  the  predictions  made  assuming  nominal  viewing 
parameters  show  very  good  agreement  with  the  available  results 
of  laboratory  experiments  in  which  neither  the  viewing  parameters 
nor  the  number  of  looks  was  constrained.  In  particular,  the 
experiments  confirm  conclusions  1,  2,  and  4 above.  Since  the 
frequency  cell  width  on  the  display  was  constant  through  the 
experiments,  they  did  not  provide  a basis  for  evaluating 
prediction  3« 
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2.0  VISUAL  PROCESSING 

The  visual  processing  systems  of  humans,  as  well  as 
those  of  other  higher-order  mammals,  are  highly  complex, 
consisting  of  millions  of  functional  parts.  These  systems 
have  been  the  subject  of  studies  spanning  decades  by  scientists 
from  several  disciplines,  such  as  psychology  and  neurophysiology. 
The  discussion  of  this  section  is  based  primarily  on  the  work 

of  neurophysiologists.  However,  the  crucial  visual  parameters 
required  for  the  performance  analysis  were  derived  from  the 
results  of  psychophysical  experiments,  in  which  subjects 
communicate  their  responses  to  physical  stimuli. 

2.1  General  DescriDtlon 

Figure  2.1  shows  the  principal  components  of  the  visual 
system.  Each  eye  includes  optical  components,  such  as  a lensr 
and  a retina  which  includes  not  only  millions  of  light  sensitive 
receptors  called  rods  and  cones  but  millions  of  other  cells 
which  process  the  outputs  of  the  receptors.  The  processed 
outputs  are  transmitted  via  neurons  comprising  the  pair  of 
optic  nerves  to  the  lateral  geniculate  for  further  processing, 
and  then  to  the  visual  cortex  for  yet  further  processing. 

The  optical  elements  of  the  eye  produce  an  image  of  the 
field  of  view  on  the  mosaic  of  photosensitive  receptors  of 
the  retina.  Even  when  the  eye  is  focused  on  a point  in  the 
field  of  view,  the  image  of  that  point  is  spread  over  an  area 
that  includes  several  receptors;  and  it  has  been  concluded  that 
the  resolving  capability  of  the  system  is  not  limited  by  the 
size  and  density  of  the  receptors.  Even  though  the  image  of 
a large  object  is  properly  focused  in  the  central  area  of  the 
retina  (the  fovea),  it  may  not  be  exactly  focused  in  other 
areas;  and  these  refractive  errors  tend  to  increase  with 
distance  from  the  fovea.  However,  the  degradation  is  gradual, 
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FIGURE  2.1  Representation  of  the  Visual  System. 
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and  for  images  of  the  size  that  will  be  considered,  the 
refractive  errors  will  be  assumed  to  be  negligible.  Further- 
more, as  will  be  discussed  later,  the  functional  units  of  the 
retina  increase  in  size  with  displacement  from  the  fovea  so 
that  the  refractive  errors  are  even  less  significant. 

The  retina  includes  a mosaic  of  photosensitive  receptors, 
the  rods  and  cones.  The  rods  function  at  low  (scotopic)  light 
intensities,  whereas  the  cones  function  at  medium  to  high 
(photopic)  intensities.  There  are  many  more  rods  (roughly 
100  million)  than  cones  (almost  10  million)  in  the  retina. 

The  density  of  cones  is  highest  in  the  center  of  the  fovea, 
and  the  density  decreases  rapidly  with  displacement  from  the 
center  of  the  fovea,  where  visual  acuity  is  greatest  under 
photopic  conditions.  When  a person  is  "looking"  at  an  object, 
at  least  part  of  the  image  of  that  object  falls  on  the  fovea. 
The  image  of  a sonar  display  viewed  at  normal  distances  will 
cover  more  than  the  foveal  area,  so  the  characteristics  of 
peripheral  vision  are  Important  to  detection  performance. 

In  dealing  with  objects  and  their  images,  it  is  convenient 
to  specify  their  extent  in  terms  of  the  visual  angles  they 
subtend.  Angular  measures  are  also  employed  to  specify  points 
on  the  retina.  The  meridian  angle  of  a point  is  the  angle 
between  a vertical  plane  passing  through  the  center  of  the 
fovea  and  the  center  of  the  lens  and  another  plane  that  passes 
through  those  points  and  the  point  on  the  retina.  The  peri- 
pheral angle  is  the  angle  between  the  line  that  passes  through 
the  center  of  the  lens  and  the  center  of  the  fovea,  and  a 
line  that  passes  through  the  center  of  the  lens  and  the  point 
on  the  retina. 
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Typically,  lofar  displays  and  bearing-time  recorders 
are  viewed  under  photopic  conditions;  therefore,  the  cones 
and  the  processing  of  their  responses  is  of  primary  interest. 
Based  on  counts  by  Osterberg,  the  density  of  cones  (number 
per  minute2)  along  the  temporal  meridian  is 

p(e)  = 5.32  exp  (-1.88e)  + 1.24  exp  (-0.150e)  (1) 

where  e is  the  peripheral  angle  in  degrees.*  The  first  term 
accounts  for  the  heightened  density  of  cones  in  the  foveal 
region.  The  contribution  of  this  term  is  practically  negligible 
above  2\  degrees. 

The  population  of  cones  can  be  divided  into  three  groups 
according  to  the  dependence  of  their  response  on  the  frequency 
of  the  incident  light.  The  response  of  each  group  peaks  at  a 
different  frequency.  The  differences  in  spectral  response  is 
the  basis  for  color  discrimination  under  photopic  conditions. 
Subsequent  discussions  in  this  report  will  not  concern  color 
discrimination. 

The  inner  portion  of  the  retina  consists  of  layers  of 
transparent  cells  that  transfer  and/or  operate  on  the  electrical 
impulses  produced  by  the  rods  and  cones.  One  important  function 
is  performed  at  the  synapses,  areas  in  which  neural  cells  are 
contiguous.  At  an  excitatory  synapse,  the  excitation  of  one 
cell  causes  an  excitation  of  another  cell,  while  at  an  inhibitory 
synapse,  the  excitation  of  one  cell  reduces  or  inhibits  the 


^Derived  from  Equation  (1)  of  D.H.  Kelly,  "Effects  of  the  Cone-Cell 
Distribution  on  Pattern  Detection  Experiments,"  Journal  oh  the.  Optical 
Society  oh  AmeAlca,  (Vol.  64,  No.  11)  Nov.  1974,  pp  1523-1525. 
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activity  in  another.  Another  important  function  is  performed 
by  the  ganglion  cells,  that  of  collecting  (or  summing)  the 
impulses  received  from  multiple  sensors.  Each  ganglion  is 
connected  to  a single  fiber  in  the  optic  nerve,  which  contains 
about  1 million  neural  fibers.  Since  the  total  number  of 
receptors  is  much  greater  than  the  number  of  fibers,  it  is 
clear  that  if  all  of  the  sensors  are  employed,  that  some 
form  of  convergence,  such  as  summation,  must  take  place. 

If  a retina  were  to  be  uniformly  illuminated,  it  would 
be  found  that  the  responses  observed  at  the  ganglions  would 
vary  with  the  rate  of  change  of  illumination.  This  implies 
that  there  would  be  no  response  to  constant  illumination, 
which  seems  contrary  to  experience.  However,  even  when  the 
gaze  is  fixed,  the  eyes  undergo  small  involuntary  rotations 
or  drifts,  which  result  in  the  retina  moving  relative  to  the 
image  of  the  field  of  view.  Thus,  a temporally  stable  image 
with  spatial  variations  of  illumination  give  rise  to  responses 
from  the  sensors  that  are  being  swept  under  the  nonuniform 
portions  of  the  image.  The  average  drift  rate  is  a few 
minutes  of  arc  per  second.  Periodically,  a very  quick  move- 
ment called  a saccade  occurs  in  the  reverse  direction  of  the 

net  drift.  The  average  size  of  a saccade  is  about  ten  minutes 
2 

of  arc.  Although  these  movements  are  larger  than  the  resolu- 
tion of  foveal  vision,  we  are  not  aware  of  them.  We  will 
return  to  this  matter  later  in  this  section. 

As  indicated  by  Figure  2.1,  the  neural  fibers  from  a 
single  eye  proceed  to  two  separate  areas  of  the  lateral 
geniculate  body,  which  in  turn  are  coupled  to  two  distinguishable 


For  further  discussion,  see  Chapter  XIV  of  Cornsweet,  Visual  Perception. 
His  Figure  14.16  is  a convincing  demonstration  of  the  insensitivity 
of  the  visual  system  to  small  rates  of  change  of  illumination. 
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areas  of  the  visual  cortex.  The  intersection  of  the  optic 
nerve  and  the  retina  contains  no  receptors,  and  hence  con- 
stitutes a blind  spot  on  the  retina. 


2.2  Receptive  Fields 

Systems  of  many  kinds  can  be  characterized  according  to 
their  responses  to  sets  of  specific  stimuli.  This  approach 
is  of  great  value  for  investigating  systems  of  great  complexity 
such  as  the  visual  system. 

By  neurophysiological  methods,  it  is  possible  to  obtain 
a great  deal  of  information  about  the  visual  systems  of 
animals  under  laboratory  conditions.  By  the  use  of  extremely 
fine  electrodes,  it  is  possible  to  monitor  the  impulse  activity 
of  neural  cells  of  anesthesized  animals  generated  by  specific 
visual  stimuli.  The  receptive  field  of  a neuron  is  defined  as 
that  portion  of  the  retina  (or  the  corresponding  visual  field) 
that  produces  a response  in  the  neuron  when  stimulated  by  light. 

Fairly  extensive  data  about  the  receptive  fields  of  cells 

in  different  parts  of  the  cat's  visual  system  have  been  obtained. 

3 

Hubei  and  Kuffler  investigated  the  receptive  fields  of  retinal 
ganglions  and  those  of  neurons  in  the  lateral  geniculate  body. 
Subsequently,  Hubei  and  Wlesel  investigated  the  characteristics 
of  the  receptive  fields  of  cells  in  the  visual  cortex. 


Recordings  were  made  from  acutely  prepared  anesthesized 
cats.  Stabilization  of  the  eyes  was  achieved  by  means  of  a 
paralyzing  drug,  and  pupils  were  dilated  with  atropine.  The 


3This  work  is  summarized  in  D.H.  Hubei,  T.N.  Wiesel,  "Receptive 
Fields,  Binocular  Interaction,  and  Functional  Architecture  in  the 
Cat's  Visual  Cortex",  Journal  oi  Physiology,  (No.  160),  1962,  pp  106-154. 
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animal  faced  a wide  screen,  and  various  patterns  of 
white  light  were  shown  on  the  screen.  Recordings  were  made 
in  the  light-adapted  state. 

The  eyes  usually  diverged  slightly,  so  that  points  on 
the  screen  at  the  center  of  vision  of  each  eye  were  not 
necessarily  superimposed,  a condition  that  allowed  mapping 
the  field  associated  with  each  eye  separately. 

In  the  retina,  the  receptive  fields  of  the  ganglions 
were  found  to  have  circular  outlines  as  indicated  by  Fig.  2.2A. 
Each  field  consisted  of  two  definable  areas:  a central 

circular  area,  and  a surrounding  concentric  ring.  Two  types 
of  fields  were  found.  In  the  "on"-center  field,  a light 
stimulus  in  the  center  area  (in  this  case  called  excitatory) 
caused  an  increase  of  the  frequency  of  response  impulses. 

In  the  concentric  outer  area,  in  this  case  termed  inhibitory, 
the  stimulus  either  reduced  the  impulse  frequency  or  terminated 
the  impulses  completely;  and  the  cessation  of  the  stimulus 
was  marked  by  a temporary  increase  of  pulse  frequency.  For 
the  "of f "-center  field,  the  positions  of  the  excitatory  and 
inhibitory  areas  were  reversed.  The  numbers  of  each  type 
of  field  observed  were  about  equal.  In  either  type  of  field, 
if  the  two  areas  are  stimulated  simultaneously,  they  react 
antagonistically  and  produce  weaker  responses  than  when  the 
areas  were  excited  separately.  The  size  of  the  fields  were 
found  to  increase  with  peripheral  angle.  Similar  fields 
were  found  for  the  cells  of  the  lateral  geniculate  body, 
which  is  further  along  the  visual  pathway. 

Hubei  and  Weisel  investigated  the  receptive  fields  in 
the  cat  cortex,  further  down  the  pathway.  With  the  exception 
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of  the  fields  of  the  incoming  fibers  from  the  lateral  genicu- 
late body,  they  found  no  circular  receptive  fields  of  the  type 
observable  in  the  retinal  ganglions  or  in  the  cells  of  the 
lateral  geniculate  body. 

The  fields  they  observed  were  classified  in  two 
categories,  simple  or  complex.  Except  for  their  shapes,  the 
simple  cortical  fields  were  similar  to  the  fields  of  the 
retina  and  the  lateral  geniculate  body  in  that  (1)  they 
embraced  distinct  excitatory  and  inhibitory  areas;  (2)  there 
was  summation  within  each  area;  (3)  responses  of  the  two  areas 
were  antagonistic;  and  (4)  the  response  to  an  arbitrarily 
shaped,  fixed  or  moving  stimulus  can  be  predicted  from  the 
configuration  of  the  excitatory  and  inhibitory  areas.  Complex 
fields  do  not  have  the  properties  listed  above. 

Figure  2.2  B-E  illustrates  the  boundaries  of  areas  of 
simple  fields  found  in  the  cat  cortex.  All  are  characterized 
by  having  straight-line  boundaries  between  excitatory  and 
inhibitory  regions.  The  orientation  of  a field  is  defined  to 
be  the  orientation  of  its  straight-line  boundaries.  The 
crosshatched  areas  could  be  either  excitatory  or  Inhibitory, 
and  either  arrangement  appears  to  be  equally  likely.  In 
Figure  2.2  A-D  the  crosshatched  areas  are  much  smaller  than 
the  remaining  areas. 

Since  for  simple  fields  there  is  summation  within  an 
area,  the  optimum  stimulus  for  fields  of  the  type  in  Figs.  B 
and  C would  be  a rectangular  bar  of  light  or  dark  for  excitatory 
or  inhibitory  centers  respectively.  The  orientation  of  the 
bar  with  respect  to  the  center  was  found  to  be  critical: 
orientation  errors  of  5 to  10°  greatly  reduce  the  response 
or  even  abolish  it.  A stimulus  of  uniform  light  over  an  entire 
field  produced  either  a very  weak  response  or  none  at  all, 
indicating  that  the  magnitudes  of  the  responses  of  the  parts 
of  the  field  were  nearly  or  actually  balanced. 
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Receptive  fields  not  having  the  properties  defining 
simple  fields  are  termed  complex.  Hubei  and  Wiesel  found 
many  complex  fields  of  diverse  character.  An  example  of 
special  interest  is  depicted  in  Fig.  2.3.  The  overall  extent 
of  this  field  is  2%°  by  3°.  The  maximum  response  to  a 
stationary  stimulus  was  obtained  by  a bar  of  light  1/8°  wide 
extending  across  the  field  in  the  direction  Indicated  by  the 
dashed  lines.  The  position  of  the  bar  was  not  critical,  but 
its  orientation  was.  As  Indicated  by  the  lower  portion  of 
Fig.  2.3,  the  response  to  a pulse  of  light  was  the  combination 
of  both  on  and  off  responses.  For  complex  fields,  such  mixed 
responses  are  not  uncommon;  this  is  in  contrast  to  simple 
fields  where  the  optimal  stimulus  provokes  either  on  responses 
or  off  responses.  The  cell  was  found  to  respond  strongly  and 
continuously  when  the  optimum  stimulus  was  moved  across  the 
field  at  a steady  rate,  and  for  this  field,  the  optimum  rate 
of  movement  was  about  one  degree  per  second.  Continuous 
response  could  also  be  maintained  by  small  to  and  fro  movements 
of  the  stimulus. 

As  discussed  previously,  the  normal  involuntary  eye  move- 
ments cause  the  retina  to  drift  slightly  with  respect  to  a 
stationary  image.  Consider  a constant  light  bar  stimulus  of  a 
specified  orientation,  size,  and  shape.  We  would  expect  to 
find  a large  population  of  fields  of  the  type  depicted  by  Fig. 
2.2B  that  at  least  approximately  match  the  specifications  of 
the  image.  As  the  retina  drifts  under  the  image,  several  of 

these  fields  will  produce  transient  responses  in  succession.  In 
contrast,  consider  a complex  field  of  the  type  discussed  above 
encompassing  the  area  of  image  excursion  whose  optimum  stimulus 
approximately  matches  the  specifications  of  the  applied  stimulus. 
As  long  as  the  stimulus  is  maintained  and  there  is  no  gross  eye 
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FIGURE  2.3  Pertaining  to  One  Type  of  Complex  Receptive  Field. 
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movement,  the  output  of  such  a cell  would  continuously  be  main- 
tained, and  would  produce  no  sensation  of  drift. 

The  type  of  complex  field  discussed  above,  if  present 
in  the  human  visual  system,  would  be  expected  to  play  an 
important  role  in  the  detection  of  elongated  signal  traces 
in  a background  of  noise  markings.  Comparable  data  on  the 
receptive  fields  of  humans  do  not  exist.  However,  there 
are  data  from  psychophysical  experiments  (discussed  in  the 
next  section)  that  can  be  employed  if  the  pertinent  hypothesis 
advanced  by  Hubei  and  Wlesel,  op.  cit.,  is  accepted. 

In  essence  that  hypothesis  states  that  a field  of  higher 
complexity  is  achieved  by  collecting  the  outputs  of  an  appropriate 
set  of  fields  of  lower  complexity. 

Two  related  examples  will  serve  as  illustrations.  A simple 
cortical  field  of  Type  B in  Figure  2.2  could  be  synthesized 
from  a set  of  overlapping  geniculate  fields  (Type  A)  whose 
inner  areas  have  diameters  equal  to  the  width  of  the  inner 
area  of  the  Type  B field,  and  whose  centers  fall  on  the 
centerline  of  the  B field. 

A set  of  these  Type  B fields  with  the  same  size  and 
orientation,  with  center  overlapping  and  deployed  with  their 
edges  along  a line  orthogonal  to  their  axes,  could  be  employed 
to  synthesize  the  complex  field  described  in  previous  paragraphs. 
To  achieve  the  mixed  temporal  response,  an  equal  distri- 
bution of  on  and  off  center  fields  would  be  required.  This 
mode  of  synthesis  raises  a question.  Suppose  the  optimum 
stimulus  were  projected  on  some  area  of  the  field?  What 
happens  if  another  bar  of  similar  size  and  orientation  but 
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displaced  from  the  first  is  projected  onto  the  field?  If 
the  separation  were  sufficient  so  that  there  were  no  inhabi- 
tory  responses  induced  in  the  pair  of  simple  fields  whose 
centers  coincided  with  the  stimulus  bars,  then  a response 
stronger  than  that  from  the  first  bar  would  be  expected. 

If  the  separation  were  reduced  so  that  the  bars  were  Just 
abutting,  the  response  would  be  less  than  maximum  because 
the  (now)  single  stimulus  does  not  have  optimum  width. 

2.3  Receptive  Field  Data 

To  predict  the  capability  of  sonar  operators  of  detect- 
ing elongated  signal  traces  in  a background  of  noise  markings, 
it  would  be  useful  to  have  the  following  information  available 

1.  The  presence  or  absence  of  complex  fields 
of  the  type  just  described  in  the  previous 
section  in  the  human  visual  system. 

2.  Given  the  existence  of  such  fields,  their 
distribution  over  the  retina;  the  overall 
size  of  these  fields;  the  dimension  and  rates 
of  movement  for  optimal  stlmulae;  and 
finally  the  response  function  over  the 
stimulus  areas. 

Due  to  the  severity  of  the  laboratory  procedures  required 
to  obtain  neurophysiological  data,  it  is  necessary  to  employ 
psychophysical  experiments  to  infer  the  data.  (Even  if  neuro- 
physiological methods  could  be  employed,  it  would  be  difficult 

to  obtain  all  of  the  data.)  With  psychophysical  experiments, 
it  is  difficult  to  deduce  the  desired  field  properties  be- 
cause the  receptive  fields  are  highly  overlapped,  and  a single 
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stimulus  may  excite  many  of  them.  The  final  decision  by  the 
observer  may  be  influenced  by  the  outputs  of  at  least  several 
fields.  Results  could  also  be  influenced  by  involuntary  eye 
movements.  However,  if  the  durations  of  the  stimuli  are  brief, 
enough,  such  effects  are  negligible. 

A search  of  the  literature  has  not  revealed  any  psycho- 
• physical  experiments  directly  pertaining  to  elongated  receptive 
fields.  However,  two  sets  of  experiments  have  been  reviewed 
that  provide  data  pertaining  to  circular  receptive  fields. 

Given  the  characteristics  of  circular  fields,  it  should  be 
possible  to  deduce  the  characteristics  of  an  elongated  field, 
given  that  the  latter  is  synthesized  from  a set  of  the  former 
(the  Hubel-Wiesel  hypothesis). 

4 

A set  of  experiments  were  conducted  first  by  Westheimer, 

c 

and  then  by  Enoch,  Sunga,  and  Bachmann  with  concentric 
circular  light  stimuli. 

The  inner  portion  of  the  stimulus  was  a small  flashing 
disc  of  light,  maintained  at  a fixed  luminance  level  and  area 
at  a predetermined  position  in  the  visual  field.  The  surrounding 
ring  was  a non-flashing  area  of  illumination  whose  level  was 
varied  to  determine  the  luminance  level  required  to  make  the 
inner,  flashing  field  appear  and  disappear.  Conducting  this 
procedure  with  different  diameters  of  the  surrounding  ring 
generated  data  pertaining  to  the  size  of  the  receptive  field 
areas . 


^G.  Westheimer,  "Spatial  Interaction  in  Human  Cone  Vision, 
J.  fhyUol,  (1967),  190,  pp.  1 39-1  54. 

5Jay  M.  Enoch,  R.N.  Sunga,  and  E.  Bachmann,  "Static  Peri- 
metric Technique  Believed  to  Test  Receptive  Field  Properties  I. 
Extension  of  Westheimer's  Experiments  on  Spatial  Interaction. 
Am.  J.  Ophthatmotogy , 70,  No.  1,  July  1970,  pp.  113-126. 
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A typical  set  of  results  from  Enoch,  et  alii,  are  shown 
by  Figure  2.4.  Here,  each  curve  pertains  to  a different  level 
of  luminance  of  the  flashing  center  disc.  A curve  is  generated 
by  selecting  the  size  of  the  surrounding  ring,  and  then  determining 
the  level  required  to  extinguish  the  appearance  of  the  flashing 
center.  The  descending  portion  of  each  curve  has  been  termed 
the  summation  arm,  and  the  ascending  portion  has  been  termed 
the  inhibition  arm.  It  is  reasonable  to  expect  that  light 
falling  in  the  inner  area  of  the  receptive  field  will  raise 
the  threshold  for  observing  the  flashing  center.  If  the  size 
of  the  surrounding  ring  is  increased,  but  not  beyond  the  summation 
area  of  the  receptive  field,  then  a lower  level  of  light  inten- 
sity will  produce  a given  amount  of  light  flux.  However,  if 
part  of  the  surrounding  ring  falls  on  the  inhibition  area  of  a field, 
then  the  threshold  for  observing  the  flashing  center  light  will 
decrease,  and  the  intensity  of  the  center  ring  would  have  to 
be  increased  to  achieve  masking  of  the  flashing  center. 

For  the  three  upper  curves  of  Figure  2.4,  the  intrusion 
of  the  surrounding  light  ring  into  the  area  of  inhibition  is  easily 
observed.  At  the  lowest  level,  the  inhibition  is  absent,  or 
poorly  developed. 

As  long  as  an  increase  in  the  area  of  the  surrounding  ring 
causes  a significant  change  in  threshold,  the  area  of  neural 
interaction  has  not  been  exceeded.  This  limiting  area  of 
neural  interaction  is  regarded  as  an  estimate  of  receptive 
field  size.  However,  as  seen  from  Figure  2.4  this  value  is 
subject  to  interpretation  since  the  rate  of  change  of  threshold 
with  ring  area  decreases  with  increasing  area. 
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Source:  Enoch,  et  alii,  op.  cit..  The  log  background  luminance 

necessary  to  just  make  the  flashing  test  field  disappear  (as- 
cending threshold)  determined  for  each  of  several  different 
size  backgrounds. 


FIGURE  2.4  Masking  Threshold  Values. 
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Data  were  obtained  with  the  flashing  spot  at  different 
periperal  angles  (the  data  of  Figure  2.4  pertain  only  to  the 
foveal  area).  Figure  2.5  shows  both  ring  areas  corresponding 
to  the  minimum  threshold  and  the  limit  of  interaction  obtained 
by  Westheimer  and  by  Enoch,  et  alii.  For  the  minimum  values, 
the  agreement  of  the  two  sets  of  data  is  good.  For  the  limit 
of  interaction,  the  data  are  somewhat  disparate,  but  this  is 
not  surprising  in  view  of  the  discussion  in  the  previous  para- 
graph, and  because  the  value  may  be  dependent  on  overall  illu- 
mination level. 

For  the  diameter  of  the  center  area  of  a field, 
Westheimer* s data  are  approximated  by  a line  given  by 

d$  = 7 + | e | / 20  , arc  minutes  (2) 


where  e is  the  peripheral  angle  in  arc  minutes.  For  the 
outer  diameter  of  the  field,  the  average  of  the  data  of 

Westheimer  and  Enoch  et  alii  is  about  three  times  that 
of  the  center  area. 

If  it  is  assumed  that  a simple  line  detector  is  synthesized 
from  a set  of  circular  fields  whose  centers  fall  on  the  same 
straight  line,  the  result  would  appear  as  shown  in  Figure  2.6 
for  a field  centered  on  the  fovea.  The  width  of  the  center 
portion  of  the  field  is  assumed  to  be 

ws  = 7 + | e | / 20  , arc  minutes  (3) 

The  width  of  the  outer  portion  of  the  field  is  assumed  to  be 
three  times  that  of  the  inner  portion. 
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Sourced  Enoch,  et  alii,  op.  cit..  Derived  dimensions  of  areas  of 
neural  interaction  (photopic)  are  presented  as  a function  of 
displacement  from  the  point  of  fixation  (ordinate).  In  addition, 
minima  falling  between  the  summation  and  inhibition  arms  of  these 
functions  are  plotted. 


FIGURE  2.5  Field  Areas. 
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FIGURE  2.6  Representation  of  a Simple  Line  Detector. 
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The  length  of  the  field  has  not  been  specified.  It  is 
assumed  that  there  are  many  fields  of  different  lengths,  and 
that  one  of  these  will  match,  or  nearly  match,  the  length  of 
the  signal  trace. 
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LOFAR  DETECTION  PERFORMANCE 


3.1  Introduction 


Lofar  equipment  is  designed  to  permit  a sonar  operator  to 
detect  the  presence  of  a narrowband  signal  in  a noise  background. 
The  equipment  consists  of  a frequency  analyzer  and  a recording 
display.  In  essence,  the  display  presents  a time  history  of 
the  frequency  distribution  of  the  input  energy.  If  the  record- 
ing medium  is  electrosensitive  paper,  the  reflectance  varies 
inversely  with  signal  intensity.  If  a cathode  ray  oscilliscope 
is  employed,  the  luminance  usually  varies  with  the  signal  inten- 
sity. A steady  narrowband  signal  appears  as  a line  of  different 
intensity  from  the  background. 


The  objective  of  this  analysis  is  to  derive  the  probability 
of  detecting  a single  narrowband  signal  by  means  of  one  of  the 
line  detector  receptive  fields.  The  elements  involved  are  shown 
in  Figure  3.1.  The  analysis  consists  of  four  major  parts: 


1.  Derivation  of  certain  statistical  measures  of 
the  Lofar  display  markings  in  terms  of  the  in- 
put voltage  statistics. 

2.  Description  of  the  retinal  image  of  the  marking 
pattern  on  the  display  surface. 

3.  Derivation  of  statistical  measures  of  the 
response  of  a receptive  field  to  the  retinal 
image . 

4 . Calculation  of  the  detection  performance  from 
the  derived  response  measures. 
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FIGURE  3.1  Elements  of  the  Problem 
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Implications  of  the  calculations  will  be  discussed  and  compared 
to  the  results  of  detection  threshold  measurements  obtained  by 
means  of  psychophysical  experiments  conducted  in  a laboratory. 

3.2  Statistics  of  Display  Markings 

Lofar  processing  has  been  implemented  with  different  circuit 
technologies,  and  the  circuit  details  do  influence  the  transfer 
characteristics.  This  subject  has  received  considerable  atten- 
tion,^ whereas  the  second  stage,  the  visual  processing,  has 
been  almost  totally  ignored.  Since  our  emphasis  is  on  the 
latter,  a reasonably  simple  representation  of  the  electronic 
circuits  is  chosen. 

The  circuit  analog  employed  to  represent  the  Lofar  equipment 
is  shown  in  Figure  3.2.  The  input,  which  could  either  be  the 
output  of  a single  hydrophone  or  a beamformer,  is  applied  simul- 
taneously to  a bank  of  narrowband  filters.  For  convenience,  it 
is  assumed  that  the  bandpass  functions  are  rectangular  and 
contiguous.  The  output  of  each  filter  is  squared  and  filtered 
to  remove  the  high  frequency  components  prior  to  sampling,  which 
is  performed  at  a rate  numerically  equal  to  the  filter  bandwidth 
(i.e.,  twice  the  half  bandwidth).  The  sampled  data  are  converted 
to  digital  form  and  are  stored  in  a memory.  These  stored  data 
are  scanned  at  a rate  required  to  refresh  the  display.  Sample 
sets  are  advanced  at  the  analog  sampling  rate,  and  the  oldest 
samples  are  discarded. 


6e.g.,  C.  N.  Pryor,  "Calculation  of  the  Minimum  Detectable 
Signal  for  Practical  Spectrum  Analyzers,"  W aval  OAdnance.  LaboAa- 
toAy  TzchnZcaZ  RzpoAt  71-92,  August  1971. 
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The  luminosity  of 
related  to  the  sampled 


a sample  marking 
voltage  by 


is  assumed  to  be 


c + kV 


i J 


(4) 


where  c and  k are  constants 

i is  the  index  for  the  sampling  time 
j is  the  index  for  the  frequency  channel 

The  first  and  second  (mixed)  moments  for  the  marking  samples 
are,  respectively 


E(Xij)  « c + k E ( V i j ) (5) 

E(X,j  Xk,>  ■ c2  ♦ c k [e  ( V ^ j ) ♦ E(Vk)  )]  ♦ k = E(V1J  »„)  (6) 

where  E represents  statistical  expectation. 

The  output  of  the  i^h  narrowband  filter  can  be  represented 
as 


U.(t)  = P.(t)  cos  2tt  f . t + Q.(t)  sin  2-rrf.t  (7) 

J J J J J 

If  Uj ( t ) is  assumed  to  be  a Gaussian  process,  then  Pj(t) 

and  Qj(t)  are  Gaussian  processes  with  the  same  autocovariance 

K(T).  If  is  squared  and  lowpass  filtered,  the  result  is 


Vj(t)  = (1/2)  [Pj(t)  + Qj (t)l 


(8) 
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The  mean  value  of  a sample  of  Vj(t)  is 

E<V1j>  ■ PNj  <9> 

where  PNj  Is  the  power  of  the  noise  In  channel  J . 

If  a signal  Is  present  in  channel  n,  then 

E<V1n>  ' PNn  + PSn  (10) 

where  Pgn  is  the  power  of  the  signal  in  channel  n 

For  the  case  of  noise  alone,  the  second  moment  function 
for  Vj  (t ) is 

E[Vj(t)Vj(t  + T)1  = PNj  + KVT)  (11) 

If  the  noise  in  the  band  has  a uniform  spectrum,  then  samples 
that  are  spaced  by  the  reciprocal  bandwidth  are  uncorrelated. 
Furthermore,  if  the  bandpasses  are  non-overlapping  but  of  the 
same  bandwidth. 


E [Vj  ( t ) V , ( t + T)]  = P2h 


(12) 


Thus,  for  the  output  samples, 

(v  (U)  t (k,l) 

E(v1j  vk1>.  " J 2 <’3) 

l 2P„ . (1,3)  * (k.D 
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The  marking  statistics  derived  from  (5),  (6),  (9)  (10),  and  (13) 
are 


E<xij> 


c + kPN  , If  = 0 


c + k(PN+  Ps)  , if  Pc,->0 


(14) 


Sj 


E<xu  xki>  * 


(c  + kPn)  , (1  ,j ) / (k,l ) 


(c  + kPN)!  t (kPN)2  , ( 1 , j ) = (k,l ) (15) 


3.3  The  Display  Image 


To  calculate  the  response  of  a receptive  field  to  the 
markings  on  a portion  of  the  display,  it  is  necessary  to  map 
the  display  image  on  that  portion  of  the  retina.  Consider  a 
region  near  a displayed  line  whose  image  would  fall  on  the 
vertical  meridian  of  the  retina.  The  most  convenient  reference 
axis  is  the  line  which  passes  through  the  center  of  the  fovea 
and  the  center  of  the  lens.  The  point  on  that  axis  whose  image 
is  in  focus  on  the  retina  is  called  the  fixation  point,  and  it 
is  assumed  that  that  point  is  the  intersection  of  the  axis  with 
the  display  surface.  The  displacement  of  an  image  point  falling 
on  the  vertical  meridian  is  given  by  the  angle  e shown  in 
Figure  3*3>  and  the  corresponding  point  on  the  display  surface 
is  a distance  l from  the  fixation  point.  The  relationship 
between  these  quantities  is  found  from  the  law  of  sines: 
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FIGURE  3.3  Viewing  Geometry. 
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sin  e " sin  (y  - e) 

where  y specifies  the  viewing  angle 

b specifies  the  viewing  distance 

Solving  for  l and  differentiating  with  respect  to  e yields, 
after  algebraic  and  trigonometric  operations  , 

dSL  _ b sin  y a A£_  (17) 

<Te  sin2  ( y - e ) Ae 

which  relates  the  two  differentials,  and  approximates  the 
relationship  between  corresponding  incremental  quantities 
such  as  the  space  occupied  by  a mark.  The  length  of  a display 
mark  in  the  time  dimension  is 

AA  = h/t  (18) 

where  h is  the  height  of  the  record  in  the  time  dimension 

t is  the  number  of  marks  in  the  time  dimension. 

The  length  of  the  mark  image  is  found  by  substituting  (18)  in 

(17)  and  solving  for  Ae . Conversion  from  radians  to  arc  minutes 

1 

gives 


, = 1 80*60  h sin2  (y  - e)  (19) 

~ rrtb  sin  y 

For  cases  in  which  the  record  is  developed  on  electrosensitive 
paper  of  indefinite  length,  the  length  of  a display  mark  in  the 
time  dimension  is 
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I 

I 


I 

I 


_ 1 


Ail  = n = s/r 


(20) 


where  n is  the  number  of  marks  per  inch 

s is  the  speed  at  which  the  paper  advances,  in/min 


r is  the  stylus  sweep  frequency,  min 


_ i 


The  width  of  a display  mark  in  the  frequency  dimension 


is 


Aw  = w/c 


(21) 


where  w is  the  width  of  the  display 

c is  the  number  of  frequency  cells 

The  corresponding  image  width  in  arc  minutes  is 


_ 1 80*60  w/c 
m it  x 


(22) 


Prom  Figure  3 . 3 and  the  law  of  sines  we  readily  obtain 


x = b sin  y 
x sin  (y  - e) 


(23) 


and  substituting  the  result  in  (22)  gives 


w = 180*60  w sin  (y  - e) 
ircb  sin  y 


m 


(24) 


Finally,  the  angular  coordinate  associated  with  the  point 
P of  Figure  3*3  can  be  derived  from  (16).  Multiplying  its  means 
and  extremes  gives 
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i.  sin  (y  - e)  s b sin  e , (25) 

and  expanding  the  left  hand  side  gives 

£(sin  y cos  e - 

Rearranging  (26)  gives 

sin  e 
cos  e 

Thus : 


cos  y sin  e)  * b sin  e 


(26) 


_ i si n y 
b + l cos  Y 


(27) 


6 * arc  tan  b ♦ I1 cos  r (28) 

The  time  period  corresponding  to  a trace  of  length  £ on  recording 
paper  is 


T = t/s  , seconds 


(29) 


Solving  for  £ and  substituting  the  result  in  (28)  gives 


„ _ +an  Ts  sin  v 
arc  tan  t,  + Ts  cos  y 


(30) 


This  relationship  gives  the  peripheral  angle  associated  with  a 
record  duration  T for  the  specified  viewing  angle  and  distance 
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The  special  case  of  an  image  centered  on  the  retina  is 
depicted  by  Figure  3.**»  and  it  is  desired  to  determine  the 
peripheral  angle.  Applying  the  law  of  sines  to  the  largest 
triangle  gives 


sin  2e  - sin  (y  - e)  (31 ) 

h y 

and  applying  the  same  law  to  the  triangle  with  sides  b and  y 
gives 

sin  v _ sin  (it  - y - e)  _ sin  (v  + e)  (32) 

y b b 

Solving  (32)  for  y and  substituting  the  result  in  (31)  gives 

sin  2e  _ sin  (y  + e)sin  (y  - e)  (33) 

h ” b sin  y 

The  solution  of  this  equation  is  found  to  be 

tan  6 ■ Fcom  [l|b2  + (h  cos  y)2  - b]  (34) 

3.4  Receptive  Field  Detection  Performance 

The  signal  detection  performance  of  a receptive  field  will 
be  predicted  by  deriving  statistical  measures  of  its  response 
to  the  display  image. 

The  type  of  field  assumed  to  be  operative  for  this  task 
is  the  complex  field  of  the  type  discussed  in  Section  2.2 
and  depicted  by  Figure  2.3.  Such  a field  could  produce  a 
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h 


h:  length  of  trace 


FIGURE  3.4  Image  Centered  on  Retina. 
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steady  response  to  a fixed  object  despite  the  drift  movement 
of  the  retina.  However,  to  simplify  the  analysis,  it  will  be 
assumed  that  the  image  is  stationary  on  the  retina,  and  that 
the  operative  field  is  of  the  simple  type  depicted  by  2.7. 

Response  Function 

If  a single  receptor  responds  in  proportion  to  the  inten- 
sity of  the  incident  light,  then  the  incremental  response  of 
the  field  to  an  incremental  illuminated  area  can  be  expressed 
as 


The  response  of  the  receptive  field  to  a set  of  markings 
can  be  written  as 


AY  = pfAAk  X 


(35) 


where  p is  the  local  density  of  cones 


f is  the  fraction  of  cones  employed  in 
the  receptive  field  times  the  sign  of 
the  response 


AA  is  the  incremental  area 


k is  a constant 


X is  the  luminance  of  the  object  corres- 
ponding to  the  area  of  the  retinal  image. 


^ (u) 


(36) 


1 «1  J 3j(u) 


37 


Report  No.  3225 


Bolt  Beranek  and  Newman  Inc. 


where  u and  v are  the  angular  coordinates 

of  a point  on  the  retina 

p(u,v)  is  the  cone  density  function 

f(u,v)  is  the  field  weighting  function 

XAj  is  the  luminance  of  mark  ij 

and  are  the  end  points  of  the  i£h. 
time  sample 

i^j(u)  and  6j(u)  are  the  end  points  of  the 

frequency  sample  at  u. 

This  formulation  assumes  that  the  field  axis  is  aligned  with 
the  time  dimension,  and  that  the  display  surface  is  uniformly 
illuminated. 

Noise  Alone  - Mean  and  Variance 

By  taking  the  expected  value  of  (36)  and  by  making  use 
of  (14)  it  is  determined  that  the  average  of  the  field  response 
to  noise  alone  is 

m YN  = E<  y) 


/w<  C ^ * ( u ) 

1 du  E I J dv  p(u,v)f(u,v)  (37) 

’«■  Vu) 


1 


0)  lP 

= k (c  + kPN)J*du  J dv  p(u,v)f(u,v) 


(38) 


a 


B 


Where  a and  w mark  the  extent  of  the  field  in  the  time  dimension 

6 and  if*  mark  the  extent  of  the  field  in  the  frequency  dimension. 
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The  double  Integral  is  the  response  to  uniform  illumination  of  unit 
intensity.  If  it  is  assumed  that  the  response  of  a uniformly 
illuminated  corttcal  receptive  field  is  zero,  then  the  in- 
tegral is  zero  and: 


m yn  = 0 (39) 

for  noise  that  is  statistically  uniform. 

The  square  of  the  response  is 


2 

Y 


• 2 

(k  ) £ 

1 


^ M u ) 

duZ  / 

j B 4 ( U ) 


dvp(u,v)f (u,v) 


x £ Jwk  dw£j**(w)  dxp(w,x)f(w,x)XijXkl  (4°) 
k ak  l 6^ ( w) 

Since  the  mean  response  to  noise  alone  is  zero,  the  variance 
of  the  response  is  equal  to  its  expected  square.  Taking  the 
expected  value  of  (39)  and  making  use  of  (15)  gives: 


°y  = U‘)  (c 


+ kPN)2£  P du£ p 

i a.  J $j(u 


^4  (U) 

J dvp(u,v)f(u,v) 


) 


1/  v C (w) 

2_j  J dwX.  J * dxp(w,x)f(w,x) 


k*iak 


39 


(equation  continued  to  next 
page) 
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+ ( k ' ) 2 j^(k  + cPN)2  + (cPn^]22  f'  du /4>J^U^dvP(u»v)f  (u»v) 

1 J Bj(u) 

X J dw  /*^W^dxp(w,x)f  (w,x) 

a1  3j(w) 


= (k1)2 


(c  + kPN) 


je.j(u) 


dup(u ,v)f (u ,v) 


■ u ^ dw  ?/  dxp(w,x)f (w,x) 
ak 


+ (k*kPN) 


dvp(u,v)f(u,v) 


(42) 


The  first  term  of  (42)  is  the  square  of  the  mean  value  which 
Is  zero;  hence. 


Oy  = ( k ' kPN ) 2 duj^j^)  dvp  ( u , v )f  (u  , v ) (43) 

1 j La i Bj(u)  J 


40 


Report  No.  3225 


Bolt  Beranek  and  Newman  Inc. 


The  double  sum  over  the  display  Image  marks  subtended 
by  the  receptive  field  would  be  tedious  to  evaluate.  The 
procedure  which  follows  will  approximate  the  double  sum  by 
an  integral.  Applying  the  mean  value  theorem  of  integral 
calculus  to  the  inner  integral  yields 


wi  r 

r du 

ai  L 

1 

(k'kp»)  & 

4»j(u)  - e j ( u ) 

p(u,vj  )f  (u.vj  ) 

(44) 


where  vj  is  a point  within  the  interval  ^j(u),  Bj(u).  Applying 
the  same  theorem  to  the  remaining  integral  yields 


2 


(45) 


(46) 


where  A “ wi  “ ai>  the  length  of  the  i52i  mark  in  the  time 
dimension;  wm(ui)  = ’J'j(u^)  - the  width  of  a mark  at  u^. 

The  inner  sum  can  be  rewritten  as: 


Where  fcmi  = - a±>  the  length  of  the  i—  mark  in  the  time  dimension; 

wm^ui^  =^j^ui^  " the  width  of  a mark  at  u^ 


The  inner  sum  can  be  rewritten  as : 
S 


= £ Vi( ui  K(ui’vj  )f  1 ui ,VJ  )}*«( ui ) 


(47) 
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which  shows  that  the  terms  are  a product  of  the  factor  in  brackets 
and  a width  increment  wmj(ui)-  This  sum  can  be  approximated  by 
one  with  a larger  number  of  terms  with  a smaller  width  increment 
of  constant  width  Av: 


' * p [wmp( U1  K(ui’vp)f1u1’vp)]4v 


(48) 


Finally,  this  expression  can  be  approximated  by  the  integral 


; = wm(  u i ) ' dv[°(  vv)f(ui'v)] 


(49) 


where  the  integration  is  carried  out  over  the  limits  of  the 

receptive  field.  If  the  axis  of  the  receptive  field  passes 

through  the  fovea,  then  the  cone  density  can  be  approximated 

» 

as  a function  of  u^,  which  for  convenience  can  be  taken  as 
the  peripheral  angle,  and  the  integration  is  over  the  width 
of  the  field  w^.: 


II  1 1 \ 

2.  .2  wf(ui) 

2i  • 1 

S * wm(  u i r 

)(u.,°)  ! 

‘ 2 wf(ui 

dyf  ( u ^ , v 1 

) 

(50) 


The  weighting  function  is  deduced  from  the  data  of  Westheimer 
and  those  of  Enoch  et  alii  (see  (3)  and  Fig.  2.7): 


f(u,v) 


( °* 

M/2, 

- (3/2)w$ (u ) 

1, 

-(l/2)ws(u) 

-1/2, 

(l/2)ws(u) 

1 0, 

(3/2)ws(u) 

V < -(3/2)ws(u) 

< v < -(l/2)ws(u) 

< v < (l/2)ws(u) 

< v <.  (3/2)ws(u) 

< v 


(51) 
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Evaluating  (50)  with  this  function  yields 

S = (3/2^  u.'  )p2(  ui  ,0  )ws(u1  ] (52) 

and  substituting  this  result  into  (46)  produces 

oY  = (3/2)  | k kPNj  I [£miwm(ul)p  ( U1 ,0  )ws(  ul) ]fcm1 


The  indicated  sum  can  be  approximated  by  one  with  uniform 
increments  smaller  than  giving 

°V  * (3/2)  ( k kPN  I [tn,qWm(UqK(Uq'°)Ws(Uq)]iU  l54> 


Finally,  this  sum  is  approximated  by  an  integral: 

<Jy  ~ (3/2)  (k  kPN  ) f dup  (u)ws(u)Am(u)wm(u)  (55) 

' 'a 


The  angular  dimensions  of  the  images  of  display  markings  were 
derived  in  Section  3*3.  Substituting  (19)  and  (24)  in  (55) 
gives 


where 


2 
° Y 


_ 3hw 
2tc 


(1 80*60k ' k 
Trb  sin  y 


(56) 


I 


d 


0J  2 

= f dup  (u)w  (u) 
a 


3 

sin 


(Y 


U) 


(57) 
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Mean  Value  of  Response  to  Noise  Plus  Signal 

The  mean  value  of  the  receptive  field  response  will  be 
evaluated  for  the  case  in  which  a signal  is  present  in  fre- 
quency bin  k whose  trace  falls  within  the  summation  region 
of  the  receptive  field.  Using  (36)  and  (14)  yields 

dvf(u.v)  (58) 


u) 

dv  f (u ,v) 

u) 


dvf (u , v ) 

(59) 

The  first  term  of  (59)  is  recognized  as  the  mean  value  of  the 
response  to  noise  alone,  which  is  zero.  The  inner  integral  of 
the  second  term  of  (58)  is  readily  evaluated  since  f(u,v)  « 1 
in  the  summation  region  of  the  receptive  field,  and  the  result 
appears  as  the  second  term  of  (59)*  The  mean  value  of  the 
response  with  signal  is  then 

0) 

my  = k'kP$  / dup(u)wm(u)  (60) 


1c  + kP»)  ? I 


mYS  * k 'l  c + kP 


ui  'Mu) 

dup(u,0)  Z /J 


i 


j 6.5  (u ) 
jjik  j 


+ k 


[c  ♦ k(pN  ♦ ps)]  f 


dup (u  ,0)  f 


M 


M 


I \ “ 'fMuj 

k ’ I c + kPM  I Z r dup(u.O)  E /J 

\ N/  i Qi  j Bj  (u ) 

k'kPs  Z Z1  dup(u  ,0)^4»j  (u)  - Bj(u)J 
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where  wm(u)  = ^j(u)  - Bj(u)  is  the  width  of  the  display  mark 
image  at  u . 


Substituting  (24)  in  (60)  gives 

180*60wk'kPs 
mYS  " Trcb  sin  y *n 


(61) 


where 

0) 

In  * / dup(u)  sin  (y  - u)  (62) 

a 

The  required  results  of  this  section  are  (6l)  for  the 
mean  value  of  the  field  output  when  the  signal  is  present, 
and  (56)  for  its  variance  with  no  signal  present. 

Detection  Performance 


In  the  last  stage  of  the  processing  sequence  depicted 
by  Figure  3*1,  it  is  assumed  that  the  output  of  a receptive 
field  is  compared  to  a threshold.  The  performance  of  this 
stage  of  processing  has  been  examined  extensively  in  the 
signal  processing  literature.  A measure  of  the  performance 
is  the  detection  parameter  d,  which  is  defined  as  the  in- 
crease in  the  mean  value  of  the  stage  input  resulting  from 
the  signal  divided  by  the  standard  deviation  of  the  input 
with  noise  alone.  For  the  case  at  hand; 

d-IHy^/Cy  (63) 

and  substituting  (6l)  and  the  square  root  of  (56)  in  (63) 
gives 
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j _ % I 2tw  *n 

The  first  factor  of  (6*1)  involves  only  the  display  parameters. 

The  integrands  of  the  pair  of  integrals  defined  by  (62)  and 
(57)  include  retinal  functions  and  the  viewing  angle,  and 
the  limits  are  determined  by  both  display  and  viewing  parameters. 
The  net  effect  of  these  parameters  will  be  assessed  by  consider- 
ing a special  case  in  the  sequel. 

For  the  case  at  hand,  the  field  output  is  a linear  com- 
bination of  a large  number  of  statistically  independent  samples, 
and  its  distribution  approaches  Gaussian  according  to  the  central 
limit  theorem.  For  such  a case,  the  false  alarm  and  detection 
probabilities  can  be  expressed  as 

PF  - N(-g)  (65) 

PD  - N(d  - g)  (66) 


where  N( 


is  the  Gaussian  probability  distribution  function 
for  zero  mean  and  unit  variance  (67) 


g 


since  mY^ 


0 


(68) 


t is  the  fixed  threshold  value. 

If  the  detection  parameter  equals  g,  then  Pp  = 0.5.  If 
this  is  regarded  as  the  threshold  of  detection,  then  the  cor- 
responding value  of  the  detection  paremeter  d,p  can  be  employed 
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in  (64)  to  determine  the  signal-to-noise  power  ratio  required 
to  achieve  the  detection  threshold: 


differential  (for  Fp  = 0.5). 

The  receptive  field  has  a property  that  is  very  important 
for  a functional  detector.  Since  the  mean  value  of  its  output 
from  noise  alone  (niy^j)  is  zero,  the  value  of  t in  (68)  required 
to  achieve  a value  of  g required  to  achieve  a false  alarm  rate 
as  determined  by  (65)  does  not  depend  on  the  level  of  the  noise. 
That  is  to  say  that  no  further  processing  (such  as  averaging 
over  a larger  area  of  the  image)  is  required  in  order  to  deter- 
mine the  value  of  the  comparison  threshold.  A detector  with 
this  property  is  said  to  be  normalized. 

The  detection  performance  [(64)  or  (69)]  is  seen  to 
depend  on  several  sets  of  parameters  and  functions  listed 
below. 

1.  Display  parameters 


(69) 


w width  of  display  (frequency  dimension) 
c number  of  frequency  cells 
h height  of  display  (time  dimension) 
t number  of  time  samples  spanned  by  h 
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2.  Viewing  parameters 

y viewing  angle,  measured  from  the  plane 
of  the  display 

b viewing  distance 

3.  Retinal  functions 

p density  of  cones  in  center  area  of 
receptive  field 

wg  width  of  center  area  of  receptive  field 

In  general,  the  dependence  of  performance  on  these 
parameters  is  fairly  complex;  however,  the  ratio  w/c  appears 
only  once,  i.e.,  under  the  radical  sign  in  (64)  or  (69). 

This  result  predicts  that  performance  improves  as  the  raster 
spacing  w/c  in  the  frequency  dimension  increases.  However, 
this  relationship  holds  only  when  the  width  of  the  signal 
trace  is  less  than  the  width  of  the  receptive  field.  The 
radical  sign  also  includes  the  raster  spacing  in  the  frequency 
dimension,  h/t.  However,  since  the  limits  of  integration 
also  depend  on  the  display  height  h,  the  significance  of 
this  spacing,  as  well  as  that  of  other  parameters  will  be 
investigated  by  means  of  a set  of  examples. 

In  the  sequel,  estimates  of  the  retinal  functions  are 
applied,  and  the  resulting  equations  involve  only  the  display 
and  viewing  parameters . 
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3 . 5 Centered  Image,  Small  Visual  Angle 

A special  case  is  the  one  in  which  the  trace  image  is 
centered  on  the  retina.  For  this  case,  the  limits  of  inte- 
gration of  (62)  and  (57)  have  the  same  magnitude,  i.e., 
a ■ -e,  and  oo  = e.  It  is  further  assumed  that  both  the 
cone  density  and  the  width  of  the  receptive  field  are  even 
functions  of  the  peripheral  angle. 

For  this  case  (62)  can  be  written  as 
e 

In  = / dup(u)  (sin  y cos  u - cos  y sin  u)  (70) 

-e 

The  first  and  second  terms  of  the  integrand  of  (70)  are  even 
and  odd  functions  of  u respectively;  thus,  with  the  assump- 
tions about  the  symmetry  of  the  functions, 

e 

In  * 2 sin  y f dup(u)  cos  u (71) 

0 

and  if  c remains  below  15  degrees, 

In  - 21  -j  sin  y (72) 

where 

I,  = / dup(u)  (73) 

0 

For  the  same  set  of  assumptions  (57)  can  be  expressed  as 
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n y ^ 


2 2 

1^  « 2 sin  Y ( I9  sin  y + I-j  cos  Y 


) 


(74) 


where 


e 2 

I2  = / dup  (u)w$(u) 
0 


(75) 


e 2 : 

Ij  = 3 / dup  (u)w  (u)u 

0 


(76) 


With  these  results,  the  threshold  signal-to-noise  ratio  is 


d - I. 


sin  y 4tw  , S 

I«  sin  y + Iq  cos  y 1 ^N 


(77) 


and  the  threshold  signal-to-noise  ratio  is 


ps 


= d. 


[^7  I2  sin2  Y + I3  cos  y 
\4tw  ^ sin  y 


(78) 


The  integrals  designated  1^,  I2»  and  were  evaluated 
using  (1)  for  the  cone  density  function  and  (3)  for  the  width 
of  the  center  area  of  the  receptive  field.  Values  of  these 
integrals  are  listed  in  Table  3.1.  For  a given  length  of  record 
or  signal  trace,  the  value  of  e,  the  limit  of  integration,  can 
be  determined  from  (34).  Corresponding  values  of  the  integrals 
can  be  determined  by  interpolating  the  tabulated  values.  In- 
spection of  the  table  shows  that  is  very  small  compared  to 
I2,  and  that  the  ratio  of  to  I2  increases  with  e.*  Further- 
more, as  seen  from  (77)  and  (78),  these  factors  are  multiplied 
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Table  3:1 


List  of 

Values 

of  Integrals 

and  Ratios 

(deg) 

0.1 

Ll 

.610 

[z 

.441 

I * I 

l3 

2 * *3 

h * 4*2 

.919 

0.2 

1 .14 

.778 

8.22-6 

1.29 

1.77 

0.4 

1.97 

1.24 

4.66-5 

0.6 

2.64 

1.55 

1.13-4 

2.12 

0.8 

3.14 

1.75 

2.01-4 

2.37 

1.0 

3.56 

1.88 

8.02-4 

2.60 

1.5 
2 . 0 

4.35 

2.10 

6.03-4 

3.00 

A A A 

4.94 

2.24 

9.84-4 

1595 

3 . 30 

2 . 5 

5.45 

2.35 

1.47-3 

3 . 56 

3.0 

5.86 

2.44 

2.11-3 

1156 

3.75 

4.0 

6.61 

2.59 

3.90-3 

639 

4.11 

5.0 

7.24 

2.73 

6.43-3 

424 

4 . 38 

6.0 

7.79 

2.85 

9.56-3 

298 

4.61 

7 0 

8.27 

2.94 

1 .33-2 

221 

4.82 

8.0 

8.68 

3.03 

1.24-2 

174 

4 . 99 

10.0 

9.33 

3.14 

2.60-2 

121 

5.27 

p j “j 

12.0 

9.82 

3.22 

3.43-2 

94 

5 . 47 

r r o 

14.0 

10.18 

5.63 

4.20-2 

78 

5 • 63 

♦Values  of  I~  are  given  in  exponential  form.  For  example, 

^ - 6 
the  first  entry  in  that  column  is  8.22  x 10"  . 
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by  cos  y and  sin  y respectively,  so  that  the  value  of  1^  could 
be  significant  for  very  small  viewing  angles  y and  large  visual 
angles  e.  It  can  be  readily  shown  that  if  calculation  errors 
less  than  0.1  dB  are  required,  that  the  value  of  can  be 
neglected  if  the  condition 

I2/I3  > 8.2/tan2  y (79) 

is  satisfied.  Table  3.2  shows  maximum  tabulated  values  of  the 
visual  angle  e which  satisfy  that  condition  for  various  values 
of  viewing  angle  y. 

Table  3:2 

Maximum  Values  of  e Satisfying  (79) 


v 

degrees 


8.2/tan2  y 


emax  tab 
degrees 


5 

62 

3 

10 

114 

6 

1 5 

264 

10 

20 

1071 

12 

The  condition  (79)  will  be  satisfied  for  most  cases  of  interest 
because  for  a given  record  length,  the  visual  angle  decreases 
with  the  viewing  angle.  Given  that  (79)  is  satisfied,  then 
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Values  of  are  listed  in  the  last  column  of  Table  1. 

At  this  point,  the  effect  of  viewing  distance  can  be 
qualitatively  evaluated.  By  differentiating  (3*0  with 
respect  to  the  viewing  distance,  it  is  verified  that  the 
tangent  of  the  viewing  angle,  and  hence  the  viewing  angle, 
varies  inversely  with  the  viewing  distance.  Moreover, 

Table  1 shows  that  varies  with  the  viewing  angle,  and 

since  (80)  shows  that  detection  performance  varies  with  this 
ratio,  it  varies  inversely  with  viewing  distance. 

3.6  Example  and  Comparison 


The  performance  of  lofars  is  often  specified  in  terms 
of  (Pg/N)^,,  where  N is  the  spectral  density  of  the  noise  at 
the  input  at  the  frequency  of  interest,  and  Pg  is  the  signal 
power  required  to  achieve  a specified  detection  probability 
for  a given  false  alarm  probability.  Since  PN  = NW,  where 
W is  the  filter  bandwidth. 


(82) 


*The  procedures  for  utilizing  (80)  or  (81)  are  discussed  in  the 
last  paragraph  of  this  report,  which  begins  on  page  56. 
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where  (Pg^N^T  is  given  either  (8l),  (78),  or  (69). 

The  variation  of  the  recognition  differential  10  log1Q  of 
(Ps/N)t  as  a function  of  the  record  length  in  time  units  is 
shown  by  Figure  3*5  for  several  values  of  the  viewing  angle 
y.  The  display  parameters  are  those  of  one  of  the  operational 
lofar  equipments,  and  the  viewing  distance  is  assumed  to  be 
18  inches.  The  curves  show  (1)  that  there  is  a declining 
rate  of  improvement  with  increasing  record  length,  and  (2) 
that  a significant  advantage  can  be  obtained  by  decreasing 
the  viewing  angle  y.  This  effect  is  utilized  by  sonar 
operators  and  has  been  verified  by  detection  threshold 
measurements . 7 

Detection  threshold  levels  of  lofar  equipment  have  been 
measured  by  means  of  psychophysical  experiments  conducted  in 
a laboratory.  Both  signals  and  noise  are  applied  to  the 
input,  and  the  signal  level  required  for  a designated  proba- 
bility of  detection  is  determined.  The  dashed  curve  of 
Figure  3*6  connects  data  obtained  in  one  series  of  experi- 
ments. In  these  experiments,  the  operators  are  not  con- 
strained to  a particular  viewing  distance  and  angle,  nor 
are  they  limited  to  a single  look  at  the  record.  Thus 
these  data  can  not  be  employed  in  a direct  comparison  with 
the  predictions  obtained  with  the  method  developed  in  the 
report.  Nonetheless,  the  prediction  for  a viewing  angle  of 
90°  at  a distance  of  18"  is  shown  on  Figure  3.6.  These 
parameter  values  are  considered  nominal  for  the  type  of 
experiment  that  was  conducted. 


^This  effect  is  also  mentioned  in  B.  J.  Pernick  and  C. 
Bartolotta,  "Signal  Enhancement  in  Acoustic-Spectral  Displays 
by  Optical -Spatial  Filtering,"  U.S.  W avy  Journal  Und.e.A.wa.£e.>L 
Acou.6ti.c4,  Vol . 18,  No.  1,  Jan.  1968,  p.  122,  footnote  4. 
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FIGURE  3.5  Effect  of  Record  Length  and  Viewing  Angle. 
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The  lofar  process  is  often  modeled  as  an  ideal  detector 
consisting  of  a filter,  squarer,  uniform  averager,  and 
threshold  detector.  The  performance  prediction  based  on 
this  analog  is  also  shown  on  Figure  3.6.  This  prediction 
is,  at  best,  4 dB  overly  optimistic,  and  the  error  increases 
with  record  length. 

Given  reasonable  assumptions  about  viewing  parameters, 
the  results  obtained  in  this  report  appear  well  suited  for 
predicting  the  performance  of  lofars  whose  raster  parameters 
have  been  specified. 

The  performance  prediction  equations,  e.g.,  (80)  show 
a factor  t/h  under  the  radical  sign.  This  quantity  is  the 
reciprocal  of  the  length  of  a display  mark  in  the  time  di- 
mension. Thus  it  appears  that  decreasing  this  length,  i.e., 
compressing  the  display  in  the  time  dimension  would  improve 
performance.  However,  the  visual  angle  e varies  with  h,  as 
can  be  determined  from  (3*0,  and  the  quantity  in- 

creases with  e,  and  it  is  not  immediately  obvious  how  per- 
formance varies  with  h.  The  results  of  calculations, 
shown  in  Figure  3.7,  indicate  that  a modest  improvement  of 
performance  can  be  achieved  by  decreasing  h.  The  improve- 
ment is  more  significant  at  higher  viewing  angles  and  longer 
record  lengths . Laboratory  experiments  with  the  same  degree 
of  compression  (25)  show  an  average  improvement  of  0.^4  dB 
for  a ten-minute  record. 

It  is  concluded  that  (80)  and  (8l)  provide  a convenient 
means  for  obtaining  the  detection  parameter  and  the  threshold 
signal-to-noise  ratio  respectively.  Table  3.1  gives  values 


♦This  model  is  discussed  in  C.  N.  Pryor,  op.  cit.,  pp.  2-13. 
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FIGURE  3.6  Comparison  of  Predictions  and  Experimental  Data. 
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FIGURE  3.7  Effect  of  Compression  in  the  Time  Dimension. 
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of  Ix  T for  various  values  of  e.  The  value  of  e for  a 
given  trace  length  and  viewing  parameters  can  be  determined 
from  (34).  To  obtain  the  probability  of  detection,  the 
detection  parameter  is  substituted  in  (66);  the  value  of  g 
required  in  (66)  is  obtained  by  solving  (65)  for  a specified 
false  alarm  probability.  If  the  threshold  ratio  is  desired 
in  terms  of  the  noise  spectral  density  (vice  its  power)  then 
(83)  should  be  employed  in  conjunction  with  (82).  Definitions 
and  units  of  the  display  and  viewing  parameters  are  listed 
below  for  convenience. 

w Width  of  record  (frequency  dimension),  inches 
c Number  of  frequency  cells  spanned  by  width  of  record 
h Length  of  record  (time  dimension),  inches 

t Number  of  time  samples  spanned  by  length  of  record 

For  paper  recorders 

h s 
t ~ c 

s Speed  at  which  paper  advances,  inches/minute 

r Stylus  sweep  frequency,  number/minute 

y Angle  between  display  surface  and  line  of  sight 
to  signal  trace  center — See  Fig.  3*3- 
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